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Abstract 

We present the results of a five-dimensional resonant amplitude analysis of the 
A+^pK~ir + system based on 946 ± 38 reconstructed decays. These data were 
produced in 500 GeV/c 7r - -N interactions by Fermilab experiment E791. We re- 
port measurements of the amplitudes for A+ decay into nonresonant pK~ir + and 
to piT°(890), A++(1232)K-, and A(1520)vr+ and we comment on other possi- 
ble resonant enhancements. This is the first complete amplitude analysis of the 
A+^pK~ir + system. We find that (54.8 ± 5.5 ± 3.5)% of the decays are nonreso- 
nant, (19.5 ± 2.6 ± 1.8)% of the decays are via the K resonance, (18.0 ±2.9 ±2.9)% 
of the decays are via the A ++ resonance, and (7.7± 1.8± 1.1)% of the decays are via 
the A(1520) resonance. We find evidence for an increasingly negative polarization of 
the A+ baryons as a function of p\, in agreement with a recent model[l] and with 
a related measurement [2] . 

Key words: Multidimensional, Resonance, A c , Polarization 
PACS: 13.30 a,13.30.Eg,13.88.±e,14.20.Lq 



1 Introduction 

Since the discovery of charm in 1974, great progress has been made in the 
understanding of charm meson decays. Charm meson lifetimes, branching 
fractions, and resonant decays have been studied extensively [3]. While charm 
baryon lifetimes and branching fractions have also been measured, these mea- 
surements are not as complete, and no information is available on the am- 
plitude analysis of any charm baryon decay system. Measurements of charm 
baryon decays yield information regarding the relative importance of specta- 
tor and exchange amplitudes. The latter are thought to play a leading role 
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in enhancing charm baryon decay rates. The most obvious components of 
the A+ — > pK~n + decays (and charge conjugate decays, which are implied 
throughout this paper) include the nonresonant pK~rr + decay as well as the 
pK (890) and A(1520)7r + two-body decays. All three of these decays can 
be described by spectator and W-exchange amplitudes. In lowest order, the 
A ++ (1232)fr~ decay can occur only via the exchange amplitude. Exchange 
amplitudes are suppressed in charm meson decays, at least at the quark level, 
because of helicity and form-factor effects. These effects are not expected to 
inhibit exchange amplitudes for charm baryons due to the three-body nature 
of the interaction. To understand fully this system of decays, as well as other 
charmed baryon decays, a complete resonant amplitude analysis is needed. 

The charm baryon and its decay products carry spin and the charm baryon 
may be polarized upon production. Previous charm pseudoscalar meson decay 
analyses have studied structure in the two-dimensional space of the decay 
product effective masses (Dalitz plot distributions), but the spin effects just 
described require five kinematic variables for a complete description. While 
this complicates the analysis, it affords greater sensitivity to the parameters 
of interest. As a by-product of the analysis, the production polarization of the 
A+, Pa c , is also measured. This analysis is the first five-dimensional amplitude 
analysis and, as such, is unique. 



2 Formalism 

We parameterize the observed decay rate as a function of the A+ polarization, 
Pa c , and of the magnitudes and relative phases of each intermediate two- 
body resonance decay amplitude. We assume that the nonresonant decay is 
described by an amplitude that is constant across phase space. The differential 
decay rate dr (or signal density S) may be expressed as 

dr ~ S(X) = (1 + PaJ ([ £ BrKK,- | 2 + | E B r(mr)a r I-if) 

r r 

+ (1 ~ PAe) (IE^KK-ur + 1 E s r (m r K,_i,_i| 2 ) (i) 

z r r 

where a r ,m,Ap is the complex decay amplitude for resonance r given m, the spin 
projection of the A c on the z-axis, and A p , the proton helicity in the A c rest 
frame. 

B r (m T ) in Equation 1 is the normalized relativistic Breit-Wigner amplitude 
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corrected for the centrifugal barrier [4]. Given the decay mode A c — > r(— > ab)c, 
%) = (-2beiki) £ 2 ^t^- r ( 2 ) 



where 



r, = r,(i)»!=£iM (3) 

g m r F r 2 (q ) 



for resonance r of angular momentum L at the reconstructed two body mass 
m T with the momentum q (and qo when m r = m ) of a daughter particle in 
the resonance's rest frame, and with resonance mass and width m and T as 
found in Ref. [3]. Using this convention, we set B r (m r ) for the nonresonant 
decay to be 1. Fl is the strong coupling factor at the appropriate decay vertex, 
and takes the Blatt-Weisskopf form as described in Table 1. Table 2 lists the 
range of the strong interaction, R x . 

Table 1 

The expressions for F used in the Breit-Wigner amplitude. 



L 


F L 





1 


1 


(i+RWr 1/2 


2 


(9 + 3RW + RW)- 1/2 



Table 2 

The values of R used in the Breit-Wigner amplitude. 



X 


Rx 




(GeV/c 2 )- 1 


iT°(890) 


3.4 [5] 


A++(1232) 


5.22 [6] 


A(1520) 


6.29 [7] 


A+ 


5.07 [8] 



In Tables 3-6, where the amplitudes « r ,m,A p (derived using the helicity formal- 
ism described in Ref. [9]) can be seen more explicitly, the direction (4> T ,6 T ) is 
the direction of the resonance, r, in the A c rest frame, using the convention 



4 



of Ref. [10]. The primed angles refer to the direction of one of the resonance's 
daughters in the resonance's rest frame. Note that the decay amplitudes for 
each resonance may have contributions to each of the four terms in Equation 
1. 




Fig. 1. Definition of angles using 7r + as an example. In both 

figures the A+ is at rest. In the first figure, which defines (9 P , 4> p ), the x-axis is along 

the direction of motion of the A+ in the lab frame and the z-axis is the polarization 

axis, normal to the plane of production. In the second figure we define <\>k-k as the 

— *o 

angle between the plane containing the K decay products and the plane containing 
the proton and the x-axis. 

Each event in the final data sample is described by five kinematic variables 
of interest (two two-body masses and the decay angles 9 P , P , and 0ktt as 
defined in Figure 1) which are determined after the pK7r reconstructed mass 
is constrained to the A c mass. We chose the quantization axis (the z-axis in 
the A c rest frame) to be normal to the A c production plane (as defined by 
Pbcam x Pa c , where Pbeam is the beam direction and p\ c is the A c production 
direction in the lab frame). The x-axis in the A c rest frame is chosen to be the 
direction of the A c in the lab frame. 



3 Experiment E791 and Data Selection 

We analyze data from Fermilab fixed-target experiment E791, which ran dur- 
ing 1991 and 1992. The data were recorded from 500 GeV/c tt~ beam inter- 
actions in five thin target foils (one platinum, four diamond) whose centers 
were separated by about 1.53 cm. The detector, described elsewhere in more 
detail[ll,12], was a large-acceptance, forward, two-magnet spectrometer. The 
key components for this study were eight planes of multiwire proportional 
chambers, and six planes of silicon microstrip detectors (SMD) before the tar- 
get for beam tracking, a 17-plane SMD system and 35 drift chamber planes 
downstream of the target for track and vertex reconstruction, and two multi- 
cell threshold Cerenkov counters for charged particle identification. 
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Table 3 

Amplitudes for the A+(± + ) -► (K*°(890)(l-) -» K-7r+)p(± + ) decay mode. 



m A p Amplitude 

\ \ E 1 e i ^i df x (fl-o ) dl (^)e< + E.e 1 ^ df_ 1 (fl-o) dJ o (^)e i0 K*° 

22 ft 2 2 ^ 

i .1 df , (0 r o) dj (^) + E 4 e^4 df _ 1 (fi-o) dL 10 (^)e i{ ^°-^) 

2 2 2 2 

_i i Sie ^i d^ , ((9-o) di (^)e- i( *r°-^ ) + E 2 e^ 2 d* ^ (0 r „) dj (^) 

2 2 2 2 

-i -\ E 3 e^ dK^o) di (^)e _i V + £46^4 d^ ^^o) dL 10 (^) e -^K 



Table 4 

Amplitudes for the A+(± + ) -► (A++(1232)(| + ) -> p^+)R- decay mode. 



m A p Amplitude 



i \ F ie ^i df ! (# A+ + )4 1 (9' p )+ F 2 e i ^ 2 df_ ! (£ A+ + ) d 1 ! x (0 p )e^A++ ~K) 

22 22 22 22 

i .1 F ie ^i df 1 (0 A++ )df_ 1 (^)e i ^ + F 2 e i ^ 2 d\_ 1 (e A++ )dl 1 _ 1 (e')e i <t>±++ 

22 22 22 22 

.1 i F 1 e l ^id^ 11 (^ A++ )df 1 (Oe-^A++ + F 2 e i ^d^_ 1 (fl A++ )dL(We- i ^ 

22 22 22 22 

. i . i Fie ^ x x (9 A++ )df_ 1 (0 p )e- i ^A++ + F 2 e i ^ 2 d^ , _ , (# A++ ) d 1 x _ 1 (OL) 

2 2 2 2 2 2 2 2 



An unrestrictive, open-charm event selection based on total transverse energy 
seen in the calorimeters was made in real time. Offline we require all A c event 
candidates to have a production and decay vertex longitudinally separated by 
at least 6 <j\, where <j\ is the error in that separation. The tracks and vertex fits 
satisfy x 2 requirements. The vertex must be formed from proton, kaon, and 
pion tracks identified as such by the Cerenkov particle identification system. 
Proton candidates are rejected when projected into regions of the detector 
with poor particle identification efficiency. To reject more background, we 
require that the magnitude of the vector sum of the transverse momenta of 
all the secondary tracks with respect to the flight path of the reconstructed 



6 



Table 5 

Amplitudes for the A+(± + ) -> (A(1520)(§~) -> pK-)ir + decay mode. 



m A p Amplitude 



H ie ^i df , (0 A (i52O)) df i (<£) + H 2 e^2 df (0 A(152O )) d 1 , , ^W-tf.) 

22 22 22 22 

-(fTie^i df x (0 A (i52O)) df _i (^)e^p + H 2 e^ 2 df _ x (# A(1520) ) d^_ x (^)e^(™)) 

22 22 22 22 

ff ie ^i din (0 A (i52O))f i (^ p )e- i ^(i5 2 o) + H 2 e^H 2 d!^! (0 A(152O) ) d 1 , , (0 p )e-^p 

22 22 22 22 

(ffie^i d^i(0 A(152 o))df_i(^)e- i ^(i52o)-^) + H 2 e^2 dK_ L (e A(1520) )dK_,(e' p )) 



Table 6 

Amplitudes for the nonresonant A+(^ + ) — > pA~7r + decay mode. 



m A p Amplitude 



2 2 H — c 

-i i JV_ + e*"-+ 



A c candidate be < 0.4 GeV/c. To eliminate reflections from D + — > K~7r + 7r + 
decays and from D + , D+ — > K~K + 7r + decays, we remove all events whose 
reconstructed K7T7T mass is within the range [1.85, 1.89] GeV/c 2 or whose 
reconstructed KKir mass is within the ranges [1.85, 1.89] GeV/c 2 or [1.95, 
1.99] GeV/c 2 . 

After this preliminary stage of the analysis, we have 998±167 A+ — > pK~7i + 
signal events and 107368±366 background events. Therefore we apply further 
requirements in order to improve this signal in an unbiased way. This is ac- 
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complished by a selection on the output of an artificial neural network. The 
network is trained using a Monte Carlo (MC) sample of A+ — > pK~7i + decays 
for signal and events from the wings of the pK~n + mass distribution in data 
for background. The variables used in the training include all those described 
in the preliminary analysis above and the transverse miss distance between 
the primary vertex and the line of flight of the reconstructed A c candidate, the 
scalar sum of the p\ of all the secondary tracks with respect to the flight path 
of the reconstructed A c candidate, the significance of the separation between 
the secondary vertex and the closest target foil edge, the calculated proper 
lifetime of the reconstructed A c , the ratio of the distance of each of the decay 
tracks from the secondary vertex to its distance from the primary vertex, and 
the minimum and product of the three ratios. 

The cut on the neural net output is chosen to maximize Ns : mc/\/ A^s,mc + A^b 
where N S:MC is the number of MC signal events (scaled to our data size) and 
Nb the number of background events in the signal region. After the cut on the 
neural net output, 2271 real events in the pK~7i + mass range of [2.18, 2.38] 
GeV/c 2 survive, as seen in Figure 2. 




2.2 2.225 2.25 2.275 2.3 2.325 2.35 2.375 

m(pKTr) (GeV/c 2 ) 



Fig. 2. A c — > pKw signal used in this analysis. 



4 Differential Decay Rate Fit 

Once the final data set is established we constrain each reconstructed pKn 
mass to the mass of the A c . This forces all candidate events, whether in the 
signal peak or in the wings, to have the same decay phase space as that of the 
signal. We describe each event by the five independent kinematic variables of 
interest (m 2 (Kir), m 2 (j97r), cos(# p ), <f> p , and (f>Kn)- 
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Next we use an extended maximum likelihood technique to fit the data. The 
likelihood is assumed to have the form 

In C = £ In A - In [2nN prcd ] + [N ^~J^" ? ) , (4) 

where N prcd = N s + N B and iV obs are the predicted and observed numbers of 
events, and £j is the likelihood of each event defined by a joint probability den- 
sity in the five-dimensional phase space of the decay and the one- dimensional 
space of unconstrained m P Kn (Fig. 2). 

The likelihood for an individual event is given by 

N s N B 
Hi = G(m pKnii ,x F )S(xi)A(xi) + Q(m pK ^ ti )B(xi) (5) 

Here G and Q represent normalized Gaussian and quadratic functions of the 
three body mass m P K-K while S(x), A(x), and B(x) are the signal, acceptance, 
and background normalized densities in the five-dimensional space of decay 
kinematic variables described previously. Note that the width of G is depen- 
dent on xp, the scaled longitudinal momentum of the A c -candidate, and is 
parameterized as a = a a x xp + &b- Our xp range is essentially [-0.1, +0.4]. The 
function S(x) is described in Equation 1, while A(x) is a model of the accep- 
tance, a five-dimensional density of the reconstructed and surviving uniformly- 
generated MC events, and B(x) is a model of the background, coming from 
a five-dimensional density of data events in the wings of the mass plot. We 
model our acceptance, A(x), and background, B(x), using the same approach. 
This common procedure is based on the i^-nearest-neighbors method[13]. The 
X-nearest-neighbors method works on the principle that the density at a point 
in phase space is proportional to N/V(r(xi)) where N is the number of events 
which occur within volume l / (r(a ? i )) which is a multidimensional sphere of 
radius r centered at Xj. This method is augmented to ensure that each of the 
one-dimensional projections of the five- dimensional density in question (A(x) 
and B(x)) matches the corresponding projection of the data set modeled (ac- 
cepted MC events and background events from the data wings, respectively) . 

Using the fitting software, MINUIT[14], we find the parameter values which 
maximize C in Equation 4. The results can be seen in Table 7. Because the 
polarization is dependent on transverse momentum[l], we split the data set 
into three equally populated regions of p\ (0.00-0.70 GeV/c, 0.70-1.24 GeV/c, 
and 1.24-5.20 GeV/c) and find the polarization for each region (Pa c ,i, Pa c ,2, 
and Pa c ,3, respectively). The polarization versus p\ is shown in Figure 3. 

Figure 4 shows the one-dimensional projections of the overall fit using the 
parameters of Table 7 and the data in the signal region. The x 2 /DOF for the 
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Table 7 

Decay amplitudes, polarization, and mass plot parameters for A c — > pi^7r from the 
MINUIT fit with statistical errors. The two parameters without errors are the only 
ones that are fixed. 



Terms 


Parameter 


Value 


pK (890) 


Ei 


0.52± 0.17 




4>E X 


-1.01± 0.48 




E 2 


0.20± 0.10 




<t>E 2 


2.35± 0.67 




E3 


0.2l± 0.10 




1 

<PE 3 


3.46± 0.42 






0.16± 0.10 




4>E 4 


5.29± 0.55 


A ++ (1232)K ~ 


Fi 


0.17± 0.07 




J. 

9F 1 


4.98± 0.41 




F 2 


0.38± 0.13 




4>F 2 


4.88± 0.40 


A*(1520)7r+ 


Hi 


0.18± 0.09 




(Pm 


5.93± 0.52 




H 2 


0.20± 0.07 




4>H 2 


-0.06± 0.55 


Nonresonant 


N++ 


0.46± 0.26 




4>N ++ 


3.48± 0.54 




N+- 


1.00 




<f>N+- 


0.00 




N-+ 


0.18± 0.15 




<f>N-+ 


0.75± 0.71 




N 


0.94± 0.45 






l.l3± 0.36 


Polarization [0 < p\ < 0.70 GeV/c 2 ] 


Pa c ,i 


0.15± 0.21 


Polarization [0.70 < p 2 T < 1.24 GeV/c 2 ] 


Pa c ,2 


-0.22± 0.25 


Polarization [1.24 < p\ < 5.20 GeV/c 2 ] 


Pa c ,3 


-0.67± 0.15 


# Signal Events 


N S 


946± 38 


# Background Events 


N B 


1324± 43 


Background Quadratic Term 




-0.98± 10.51 


Background Linear Term 


k 


l.34± 0.48 


A c Mass (GeV/c 2 ) 


m 


2.29± 0.00 


A c Width (MeV/c 2 ) 


cr a 


20.l± 4.8 






9.3± 0.6 
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Q. 

c ~0.8 - 

o 

N 

O0.4 - 

o 

^0.2 - 

- 

-0.2 - 

-0.4 - 

-0.6 - 

-0.8 - 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 

P, (GeV/c) 

Fig. 3. The polarization of the A c as a function of the A c 's transverse momentum. 
The vertical bars represent the error as found by MINUIT. They are placed at the 
average p\ value for that region. Bin sizes are as in Table 7. The dotted line is the 
value of the polarization integrated over pj,. 

agreement of the one- dimensional projections of the model and the data is 1.06. 
Although this agreement is reasonable, there are discrepancies, especially in 
the pK-mass-squared distribution. A model with a spin | resonance decaying 
to pK and having a mass of 1.556±0.019 GeV/c 2 and width of 279±74 MeV/c 2 
reduces the overall x 2 °f the fit. However, there is no known resonance with 
this description. Trying another fit model, we find some evidence for an upper 
tail of the A(1405)(— > pK) [15]. Similarly, an improvement may be obtained 
by inclusion of nonresonant contributions beyond s-wave. In the case of each of 
these other models, the data are not sufficient to provide compelling evidence 
and/or to give significant measurements of the required additional parameters 
for the models. Figure 5 shows the one-dimensional projections of the model's 
resonant and nonresonant components, background, and total fit. 

We calculate the fit fraction associated with each decay using Equation 6. 
A fit fraction F{ is the fraction of our signal described by decay via the i th 
contribution and is the ratio of the magnitude squared of the amplitude term 
describing this contribution and the complete probability of the decay: 



J Em,A 


p 




2 


dx 


/ Em,A p 


Ej B j( m jKm,A 


p 


2 dx 



(6) 



using the notation of Equation 1 where / dx implies integrating over the decay 
phase space. 

The major systematic errors on the fit fractions come from uncertainties in 
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Fig. 4. Projections of the data within the m V K-K range of (2.265 - 2.315 GeV/c 2 ) 
and overall fit superimposed. The data are shown as a solid histogram and the fits 
as dashed lines. There are 50 bins in each plot. 



the performance of our Cerenkov counters and drift chambers, in the MC 
model for the production of A+ baryons, and in the adjustment of the K- 
nearest-neighbors model to the one-dimensional projections. In each case, a 
different model of acceptance is studied and the change in the central values 
of parameters is taken to be an estimate of the systematic error. Table 8 lists 
the uncertainties in the fit fractions due to these various sources. 

The fit fractions for A+ — * pK~ir + are listed in Table 9, along with the statis- 
tical (first) and systematic (second) uncertainties. A check of these results is 
made by fitting the data assuming a single, average polarization. The resulting 
polarization is Pa c = -0.09±0.14 (shown as the dotted line in Figure 3). The 
component fit fractions do not change significantly. 
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A ++ (1232)K" 
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3 



mass 2 (pK) (GeV/c 2 ) 2 



••• 



Background 

A(1520)ji + 



40 

20 



• ••. 



1.5 2 2.5 3 

mass 2 (p7r) (GeV/c 2 ) 2 




•«, • 



-U.D 



cosftl 



?n •/.•• • •••••• 



20 
10 



Fig. 5. Projections of the overall fit components scaled to match the number of data 
events within the m.{pKir) range of (2.265 - 2.315 GeV/c 2 ). There are 50 bins in 
each plot. 
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Table 8 



The systematic errors on fit fractions. 





Drift 


Cerenkov Acceptance Production Combined 


Mode 


Chamber (%) Counter (%) Adjust (%) Model (%) Error (%) 


pK*° (890) 


1 1 
1. ( 


U.D U.o U.U l.o 


A [LZ6Z )J\ 


9 ^ 


19 n 7 fl 1 9 Q 
i .z u.i u. i z.y 






n ^ n a n 1 i i 

U.O U.^ U.l 1.1 


INFoTI VP^OTl P TIT 


3 9 


11 HQ D 1 ^ 

1.1 U.a U.l 0.(J 


rp i i /■ "\ 

Table 9 

The fit fractions 
final fit. 


for A+ — > pK~ir + 


with statistical and systematic errors from the 


Mode 


Fit Fraction (%) 




pK*° (890) 


19.5±2.6±1.8 




A++(1232)K- 


18.0±2.9±2.9 




A(1520)vr+ 


7.7±1.8±1.1 




Nonresonant 


54.8±5.5±3.5 




Our results are compared to other published results in Table 10. The compo- 
nents in this decay do not demonstrate significant interference. This is most 
easily demonstrated by the fact that the fit fractions in Table 9 sum to near 
unity. Although we have better statistics than previous measurements, the 
uncertainties are comparable due to our more general fit. 


Table 10 

A c branching ratios relative to the inclusive A+ — > pK~ir + branching fraction. The 
NA32 and ISR values were calculated from one-dimensional projections only. 


Mode 


E791 


NA32[16] ISR[17] 


piT° (890) 


0.29±0.04±0.03 


0.35±g;g?±0.03 0.42±0.24 


A++(1232)K- 


0.18±0.03±0.03 


0.12;°;°^±0.05 0.40±0.17 


A(1520)vr 


0.15±0.04±0.02 


0.09lg;^±0.02 


Nonresonant 


0.55±0.06±0.04 


0.56±g;g|±0.05 



5 Discussion and Conclusions 

Significant resonant and non-resonant branching fractions are found in this 
analysis of 946 ± 38 A+ — > pK~ir + signal events, the largest sample ever 
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analyzed. The size of the sample allows for inclusion of relative phases of the 
various contributions and a full accounting of spin and production polarization 
for the first time in such an analysis. The A(1232) ++ A'~ and A(1520)7r + decay 
modes are seen as statistically significant contributions for the first time, even 
when uncertainties associated with phases and other variables are included. 
The observation of a substantial A ++ K~ component provides strong evidence 
for the W-exchange amplitude in charm baryon decays. The observed com- 
ponents of the A+ — ► pK~7i + decay do not interfere significantly. Finally, we 
find no evidence for either A(1600)7r + or S(1660)7r + in the data. 

In spite of the improvements in the A+ — > pK~ir + results coming from in- 
clusion of the newly observed contributions, there remains poor agreement 
between the fit and the data in the pK-mass-squared projection. Additional 
data from new experiments are needed in order to conclusively demonstrate 
additional resonances (or their tails). 

We find evidence for an increasingly negative polarization of the A c baryons as 
a function of p^, in agreement with a recent model [1] and, for large negative 
decay asymmetry, the measurement of Ref. [2]. As suggested by the authors 
of the model, the assumption of T-invariance leads to further constraints on 
resonant phases. Again, our data are consistent with these constraints but 
additional data from new experiments are needed for a stringent test of CP 
violation. 
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